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Momentum transfer due to Coulomb interaction between two parallel, two-dimensional, narrow, 
and spatially separated layers, when a current /drive is driven through one layer, is studied in the 
presence of a perpendicular magnetic field B. The current induced in the drag layer, /drag, is evalu- 
ated self-consistently with /drive as a parameter, /drag can be positive or negative depending on the 
value of the filling factor v of the highest occupied bulk Landau level (LL). For a fully occupied LL, 
/drag is negative, (i.e., it flows opposite to /drive), whereas it is positive for a half-filled LL. When 
the circuit is opened in the drag layer, a voltage AVdrag develops in it; it is negative for a half-filled 
LL and positive for a fully occupied LL. This positive AVdrag, expressing a negative Coulomb drag, 
results from energetically favored near-edge inter-LL transitions that occur when the highest occu- 
, pied bulk LL and the LL just above it become degenerate. 

(]} • PACS numbers: 73.20.Dx 

I. INTRODUCTION 

Recently the transresistance Rt between two parallel two-dimensional (2D) systems [Q-Q has been studied exten- 
sively due to advances in measuring || techniques. In most previous studies Q-|5|] only wide systems were considered, 
in which edge effects can be ignored. However, edge effects become prominent in narrow systems especially when 
a magnetic field B is present and the highest bulk LL is completely occupied. The electronic edge states are very 
different from the bulk states. For instance, if a weak random impurity potential is present, the bulk states are 
disordered R] and occur in a series of energy bands of finite width T n (k), centered about the LLs with no states in 
the region between them. On the other hand, the edge states are degenerate with respect to the LLs. This can lead 
to profound differences in transport properties between the case when the highest bulk LL is completely occupied and 
that when it is not. 

In this paper we study the influence of edge states on the drag by solving self-consistently the Coulomb-coupled 
Schrodinger equations for two narrow Hall bars. This goes substantially beyond the simple, not self-consistent treat- 
ments of Ref. H that used a parabolic confining potential. Here the potential is evaluated self-consistently from an 
initial square-well potential. We find that a current Ii is always induced in the direction of /drive and the Coulomb 
interaction lifts the degeneracy of all occupied edge LLs, but not that of unoccupied and neighbouring-occupied edge 
LLs. The Coulomb drag obtained is two-to-three orders of magnitude larger than that in zero field due to the increase 
in the available density of states when an external magnetic field is present. As the filling factor v of the highest LL 
approaches 1, near-edge inter-LL transitions occur between this LL and that just above it. These transitions make 
^drag or AVdrag change sign when the circuit in the drag layer is closed or open, respectively. Thus, this negative drag 
£h [ is neither due to a thermal gradient S nor to to an electron-hole coupling. 

In the next section we present the formalism. In Sec. Ill we present and discuss the numerical results. We conclude 
' with remarks in Sec. IV . 



II. FORMALISM 

A. Coupled Schrodinger equations 

Consider two quantum Hall bars parallel to the (x, y) plane, separated by a distance d along the z axis, of thickness 
zero, width L x = w, and length L y = L 3> l c . The electrons are confined along x by infinitely high potential barriers. 
In a field B = — B z, the electron wavefunction in the drive layer a, when the circuit in the y direction is closed, has 
the form ip^ k (x) e tky j\[L and ip a obeys 
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The corresponding wavcfunction ^ k (x) in the drag layer b obeys 



+ (k- x/tiY - e*'4>ab(xm k (x) = E a nk r nk (x). (1) 



+ ( k - e* 2 <j> ba (x)}€ k (x) = KMx). (2) 



2m* dx 2 

Here n is the LL number, m* the effective mass, £ c — yJh/eB the magnetic length, e* = ^p-e/y^e, and e the dielectric 
constant. For simplicity spin is neglected. The Coulomb potential a & is given by 



<t> ab (x) = -2 / dx' [6p a (x') In \x - x'\ + 5p b (x') In V(z - x') 2 + eP], (3) 



As for the charge densities p a and p b , they are given by 



in the drive layer M and by 



A»0«0 - ^ E / ° dfc 5n(fc)IC fe (^)| 2 (4) 
Z7r n J - fe o 

P^ X ) = ^Y,J k dk 9 b n (kM b nk (x)\ 2 f(E b nk ). (5) 



in the drag layer. In Eqs. (3)-(5) 8p a {x) = p a {x) — po a , a — a or b, po a is the background charge density, / is the 
Fermi function, and ko = w/2£ 2 . Notice that Sk = in Eq. (^|) since no current flows through in the drag layer. We 
take the effective background densities equal poa = Pob and constant. The Schrodinger equations for the drive and 
drag layers are solved self-consistently. 

The weight function g®(k) depends on E® k and expresses the degeneracy of the LLs. When the n-th LL is completely 
filled, we have g%(k) = 1. When it is partially filled, we determine g^(fc) self-consistently Jl^) from E£ k and the level 
broadening r^(fc) due to scatterers from the assumptions (i) r°(fc) is independent of k, and (ii) the total neutrality of 
the Hall bar and the local neutrality at its center are maintained. The meaning of these assumptions becomes clear 
if we consider two limiting cases for the average filling factor of the highest LL v ^ 1. First, for an infinitely wide bar 
without edges, the energy levels E^ k is definitely degenerate with respect to the wave vector k and everywhere we 
have the same filling factor v = v. Secondly, for a narrow Hall bar without the k degeneracy in E" k and g"(k) — 1, all 
k states are occupied if v = 1. But if v < 1, not all k states are occupied and electrons prefer those low- lying energy 
states near the center of the bar. Therefore the actual width of the Hall bar, defined by the density of electrons, will 
shrink to roughly VW and in it the filling factor is still 1 even though v is less than 1. In reality, there is scattering 
which broadens E" k by, say, T"(k). Within T"(fc), we assume that each k has the same electron occupancy, so the 
local filling factor depends on the local density of states at k. This average electron occupancy depends on the density 
of states N n (k) and is reflected in 

= ^^TtS t 1 " ©(^(O) - N n (k))] 
N n {k) 

+e(P n N n (0) - N n (k)) (6) 

where Q(x) is the Heaviside step function and v n is the filling factor of the nth LL. Consequently, the electron density 
is nonzero within the width of the bar with local filling factor v ^ 1. At the center of the bar, which corresponds 
roughly to k = 0, we have gnik) — v n . This ensures the neutrality at the center which is in contrast with the 
single-layer treatment of Ref. Q where only the total neutrality was preserved and only the v = 1 situation could be 
handled. In our case, this local neutrality at the center ensures that E nk = (n + l/2)huj c at k = 0. In Eq. (Q), Sk || y 
is the shift of the Fermi surface || that results from the application of a current /drive || y through the drive layer in 
the presence of scattering. 



B. Currents 



The drive and induced currents, /drive and Ii, are given by 
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/drive = -low 2 J2 I dx I . dk 9n(k)(k - x/%) |< fc (x) \ 2 f(E^ k+sk ) , (7) 



and 

h = - W I dx f° dk g b n (k) (k - x/il)\^ b nk {x)\ 2 f{E b nk ). (8) 

where Jo = eh/2irm*w 2 . Ii consists not only of the current induced by momentum transfer but also of the classical 
E x B-drift [[Ll] current in the direction of Bx < E > where < E > is the average of a finite electrostatic field E 
exerted on the drag electrons by the charges which accumulate at the edges of the drive layer and produce the Hall 
voltage. This field E opens up the circular orbit of the drag electrons by accelerating those moving opposite to it and 
decelerating those moving in its direction jll| . Thus it induces an average current I es in the direction of /drive parallel 
to Bx < E > given by 

I es = -I w 2 [dx f° dk g° n (k) (k - x/£ 2 c )\ X nk(x)\ 2 f(E° nk ), (9) 

and Xnk{x) obeys Eq. (1) with a — > b, <f> — > <fi, p — > p, and E — > E. 4>{x) is given by Eq. (2) with the changes p a — > pj,, 
Pb — > p a , and 

P>(x) = / ° dk 9 a n(k)\Xn k (x)\ 2 f(E° nk ). (10) 



The use of the equilibrium energy E® k in the drag layer to derive <?°(fc) and of the Fermi function in Eq. @ eliminates 
electron transitions between different fc's. Thus the current due to momentum transfer is 

/drag = Ii — Ies- (11) 

For fully occupied LLs we have g^(k) — g^(k) — 1 and /drag vanishes but the classical < E > x B-drift I cs does not. 

C. Explicit relation between the drag current and drag voltage. 

When the circuit in the drag layer is opened, apart from the voltage induced by the < E > x B-drift, a drag voltage 
AVdrag = LId iag At8(z) I eW develops along the bar as shown in Fig. 1. In a short time At it produces an intermediate 
current, 

/?0 e AVd rag AT/ ("max + v)(? n 

/i„t = g = AT/ drag — (12) 

across the width, due to the Edrag x B-drift. This results in an electric field Z?i n t = Ii n t8(z)St/eL, across the width, 
which produces another Ej nt x B-drift current opposite to /drag- The steady-state is reached when these two currents 
are equal so that all the E x B-drifts are balanced. Thus, E- mt is given by 

p _ e 2 (n max + t/) 2-kK 

Eint - r Al/ drag At 6{Z) - -z— ; — T-'drag (loj 

with Umax the highest occcupied bulk LL number. The drag voltage is then 

2-Kfi L 

AVdrag --7 ; /drag, (14) 

("-max + v)e A W 

where the minus sign is inserted according to the usual Ohm's law convention; /drag is given by Eq. (9). 
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III. RESULTS AND DISCUSSION 



We consider a GaAs sample at zero temperature, with d — 5£ c , w — 20£ c , and B chosen such that £ c = 100A. The 
Schrodinger equations for the drive and drag layers are solved self-consistently and the average electron and current 
densities of the drag-layer, with n max = 1 and v = 0.5, are plotted in Fig. 2. The solid and dotted curves represent 
equilibrium and nonequilibrium quantities, respectively. Fig. 3 shows that the current density of the main figure is 
not symmetric. E® k and <7°(fc) are shown in Fig. 4. The rapid change of g"(fc) near the edges reflects the structure 
of the boundaries between edge and disordered bulk states. As shown in Fig. 2, the current-density (jb) plot has 
qualitatively five different zones, namely the outer-edge (I0£ c > \x\ > 6£ c ) and inner-edge (6£ c > \x\ > 3.5£ c ) zones 
on both sides of the bar, and one bulk zone (3.5£ c > \x\ > 0). In the outer-edge zone jb is opposite to that in the 
inner-edge zone and in the bulk zone it is negligible. 

To understand the results of Fig. 2 we consider the single-electron motion. We first treat the equilibrium case 
(/drive = 0, 5k = 0). For finite B, jb near the edge is affected by the sharp changes of the electron density n(x) and 
of the effective confining potential. Near the left edge, the gradient Vn(x) is so large that there are more electrons, 
with guiding center at x + Sx, e.g., near the maximum of n(x) at x < —9£ c , flowing into the paper (y direction) than 
adjacent electrons, with guiding center at x — Sx, flowing out of the paper; this gives rise to a local jb, at x, in the 
negative y direction. This diamagnetic-drift current is thus proportional |lf[] to Vn(x). In the same region the local 
effective electric field E c ff = — W c r, produced by electrons in the center of the bar, pushes the electrons towards 
the edge. This produces a E ff x B-drift current ]fl} | pointing in the same direction as the diamagnetic-drift current. 
However, for electrons located further away from the edges E e ff changes direction. It pushes the electrons away from 
the edge towards the center of the bar, and produces a E c ff x B-drift jb, proportional (ffj to W e g(x) — — E c g that 
points in the positive (negative) y direction at the left (right) edge. In most of the outer-edge region, the E e ff x B-drift 
is comparable to the diamagnetic-drift whereas in the inner-edge region the E c ff x B-drift dominates. As for the 
bulk zone, Vn(x) and W c a{x) are almost zero and so is jb- If we define the average drift velocities carried by the 
outer-left-edge current, inner-left-edge current, inner-right-edge and outer-right-edge currents as —v l out , v l in , —v{ n and 
v^ ut , respectively, with corresponding average electron densities n l out , n\ n n\ n and n r ont , the total current density is 

h ~ «-out u out — n out v out + n \n v \n ~ n ln v in (15) 

When no current is driven (/drive = 0), we have n r out = n l out , n T in = n| n , i£ ut = v l out , v~[ n = v\ n , or in other words 
the current density on the left edge balances out that on the right edge and thus h = 0. When /drive is switched on 
(5k ^ 0), E h nk and g h n (k) change. This affects the drag electrons in two ways. Firstly, on the average there are more 
electrons at x > than at x < due to <j) b (x), i.e. n r m + n r ont > n[ n + n l out and this leads to an imbalance of jb- 
Secondly, the inner-edge zone, shown in Fig. 2, expands into the bulk zone due to the change in g^(fc). The resulting 
total /drag is parallel to /drive for v = 0.5. As v increases, the inner-edge zone expands further, due to the increase of 
g^(k) in the bulk zone (the guiding-center coordinate is proportional to k), as contrasting Fig. 2 with Fig. 5 shows; 
There is an increase in both n\ n and n[ n with n- n < n\ n and thus decreases /drag- for u integer. 

As v approaches an integer, g^(k) — > I in the bulk is the same as that in the edge zones and there is no bulk activity 
contributing to the drag current. However, the eigenvalues E h nk , shown in Fig. 7, with n = n max near kl c = 9 and 
n = n mH + 1 near kl c = 7, are degenerate. Now the electrons prefer to occupy the inner-edge zone, lower-energy sites, 
with n max + 1 and 4.8 < \k2£ c \ < 7.6, rather than the outer-edge zone, higher-energy sites, with n max and \ki£ c \ > 8.8. 
This results in the current density shown in Fig. 6. Since the drag current is the total induced current minus the 
classical < E > x B-drift, we view these empty states as holes in the n max th LL responsible for /d rag - These holes 
are responsible for the outer-edge currents. Since there are more holes at the left edge than at the right edge and 
holes are moving opposite to the outer-edge electrons of the drive-layer, the induced current flows opposite to that 
due to electrons. On the other hand, electrons in the inner-edge zones of the n max + I LL have a total current in the 
same direction as that due to the holes.. When the sum of these two currents exceeds that due to electrons in the 
bulk n max th LL, the total current /drag, calculated by integrating jb(x) over the Hall bar, flows opposite to /drive in 
contrast with the B = case where it flows in the direction of /drive for two electron layers. Accordingly, the drag 
voltage, shown on the right axis of Fig. 5, changes sign when v approaches I and has a minimum in the vicinity of 
v = 1. This change of sign occurs also in /drag, when allowed to flow, as Eq. (12) shows. 

The predicted negative Coulomb drag can best be tested in an open-circuit configuration, i.e., by measuring the 
voltage difference along the bar to avoid the effect of scattering in the drag layer. A transient measurement should 
be conducted to resolve the classical < E > x B-drift induced voltage from the momentum-transfer induced AVdrag 
since the response time of the classical drift is far shorter than that of AVdrag- An overshoot of the measured voltage 
would signal the negative Coulomb drag. 

Finally, we notice that the induced drag, for v integer or half-integer, is about two-to-three orders of magnitude 
larger than that at zero field B. This is due to the fact the B — states condense into LLs when B ^ and agrees 
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with the results of Ref . [Q . 



IV. CONCLUDING REMARKS 

We have shown, within a self-consistent Hartree approximation, that the current density in Coulomb-coupled narrow 
Hall bars has different current zones that change with filling factor v. In addition to the zero-field momentum-transfer 
current /drag, the Hall voltage developed in the drive layer gives rise to a classical < E > xB-drift current / os when 
the circuit is closed. As v increases from 0.5 towards I, /drag or AVdrag decreases and, when v — > I, it changes sign. 
This change occurs because electrons make energetically favorable near-edge inter-LL transitions. We expect that the 
values of v where the drag changes sign depend only weakly on the initial confining potential. In the present study 
we took the latter square-well in form and evaluated the resulting potential self-consistently. Another choice would 
be an initially parabolic potential, but the results would be qualitatively the same. 

It is evident that the above negative drag is neither due to a thermal gradient || nor due to the conventional 
electron-hole coupling since we are dealing with electron layers. Also, since we have neglected tunneling between the 
layers, it cannot be identified with the observed negative drag of Ref. As we suggested in Sec. IV, it could be 
tested with time-dependent measurements. 

Finally, it should be noted that since we considered an applied current in the drive layer in the presence of scattering, 
i.e., a clearly nonequilibrium case, the resulting drag is a dissipative one and not the non-dissipative equilibrium drag of 
Ref. H . It is a result of the broadening of the Landau levels introduced by scattering coupled with the nonequilibrium 
effect as is evident from Eq. (||) in which the occupation is determined self-consistently. The scattering is embodied in 
g n (k). If we were to treat this system as a non-dissipative one, g n {k) would be 1 and thus /drag = 0. As for the use of 
5k in the Fermi function, to account for the nonequilibrium distribution function, it should be said that it was made, 
in the spirit of Ref. ||, in order to simplify the heavily involved self-consistent calculations. A thorough treatment 
would have to determine the distribution function from the Boltzmann equation. 
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FIG. 1. Schematics of a quantum Hall bar, of length L and of width W, in a perpendicular magnetic field B. The currents 
/drag and lint are explained in the text 

FIG. 2. Equilibrium (solid curves) and nonequilibrium (dotted curves) densities and current densities as a function of x for 
rimax = 1, v = 0.5 with jo = eLu c /2ir£ c and p a = 1/2tt£ c . The dashed curve is the effective potential V e ff • 

FIG. 3. Nonquilibrium current density as a function of x for n max = 1, v = 0.5. The solid curve refers to the left-bottom 
scales and the dashed curve to the right-top scales. 

FIG. 4. Equilibrium LL energies and the weight function as a function of k for n m ax = 1, v = 0.5. The dashed line is the 
chemical potential. The dots correspond to g1(k) and the crosses to g\(k). 

FIG. 5. Equilibrium (solid curves) and nonequilibrium (dotted curves) densities and current densities, as a function of x, 
and nonequilibrium weight function (crosses) as a function of k for n max = 1, v = 0.8. 

FIG. 6. Left-bottom scales: Equilibrium (solid curves) and nonequilibrium (dotted curves) current densities as a function of 
x for 1. Right-top scales: The drag voltage AVdrag (crosses) as a function of v with n ma x = 1 and /drive — 0.08/M. 

FIG. 7. Nonequilibrium LL energies as a function of k for 72 max — 1, v — 1. The dashed line is the chemical potential. 
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